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Abstract 


The  effects  of  atmospheric  turbulence  on  millimeter  wave  propagation  are  not  as  well  understood  as  the 
corresponding  effects  on  optical  propagation,  generally  because  of  the  strong  dependence  of  turbulence 
effects  on  the  absolute  humidity  structure  parameter  Cq ’  (as  opposed  to  just  the  temperature  structure 
parameter  CT~  and  the  cross-correlation  C-jq)  in  this  frequency  range.  Scattered  results  at  35,  94,  14u,  and 
220  GHz  are  available,  but  in  almost  all  cases,  available  atmospheric  data  are  Inadequate,  generally  because 
turbulence  measurements  were  obtained  incidental  to  other  propagation  experiments.  This  paper  attempts  to 
compare  available  results  to  theory,  and  shows  that  agreement  in  most  cases  is  plausible.  An  experiment 
designed  to  characterize  millimeter  wave  turbulence  at  several  frequencies  of  Interest,  while  at  the  same 
time  determining  values  of  appropriate  atmospheric  parameters,  will  be  discussed.  Included  in  tne  planned 
investigation  are  measurements  of  the  mutual  coherence  function  showing  angle-of-arrival  effects  and 
intensity  fluctuations. 


Introduction 


The  recent  Increase  in  interest  in  the  use  of  millimeter  waves  for  both  civilian  and  military  applications 
has  caused  emphasis  to  be  placed  on  the  study  of  the  effects  of  the  atmosphere  on  radiation  in  this  frequency 
range.  Attenuation  by  oxygen  and  especially  water  vapor,  scattering  by  aerosols,  and  perturbations  by 
atmospheric  turbulence  are  among  the  subjects  being  studied  at  several  laboratories  in  this  country  and 
throughout  the  world.  One  would  expect  all  of  these  effects  to  be  different  for  millimeter  waves  than  for 
either  the  microwave  or  the  infrared/visible  spectrum  because  of  the  peculiar  problems  associated  with  the 
propagation  of  radiation  of  Intermediate  wavelengths  through  the  atmosphere. 


In  considering  atmospheric  turbulence  at  millimeter  wavelengths,  for  example,  the  contribution  of  the 
absolute  humidity  fluctuations,  measured  by  the  structure  parameter  Cq-,  must  be  considered  in  addition  to 
the  fluctuations  in  temperature,  characterized  by  the  structure  parameter  Ct~.  Furthermore,  tne  cross- 
correlation  of  these  quantities,  C^,  which  may  be  either  positive  or  negative,  must  also  be  included  in  the 
analysis  of  millimeter  wave  turbulence  effects.  These  effects  Include  fluctuations  of  both  Intensity  and 
angle-of-arrival,  and  both  of  these  phenomena  will  be  treated  in  some  detail  in  this  paper. 

In  the  fall  of  1978,  millimeter  wave  propagation  measurements  were  made  by  Georgia  Tech  at  both  94  and  IvG 
GHz  at  White  Sands  Missile  Range,  New  Mexico,  over  a  2  km  path1.  Humidity,  temperature,  and  the  visible 
wavelength  index  of  refraction  structure  parameter  Cn-  were  also  measured.  Significant  intensity 
fluctuations  were  observed  at  both  frequencies  during  these  tests.  These  fluctuations  could  be  explained 
plausibly  by  a  theory  of  turbulence  effects  developed  by  Armand,  et  al.2,  who  Included  the  effects  of  the 
structure  functions  of  both  temperature  and  humidity  in  their  calculations,  although  the  cross-correlation  of 
these  parameters  was  not  adequately  treated.  A  treatment  including  this  cross-correlation  has  been  given  by 
Hill,  et  al.  ,  whose  theory  also  gives  plausible  agreement  with  the  White  bands  results,  as  will  be  shown 
later.  Subsequent  measurements  were  made  during  winter  propagation  tests  near  Burlington,  Vermont  at  fre¬ 
quencies  of  35,  95,  140,  and  220  GHz^1’^.  Again,  little  supporting  turbulence  related  meteorological  data  are 
available  because  turbulence  results  were  obtained  incidental  to  the  propagation  measurements.  However,  the 
observed  effects  of  turbulence  showed  good  agreement  with  the  expected  theoretical  dependence  of  Che  log 
amplitude  variance  on  both  frequency  and  range,  as  will  be  discussed  later.  The  results  of  these  measure¬ 
ments,  together  with  a  discussion  of  their  comparison  to  the  theory  of  Hill,  et  al.  is  the  main  thrust  of 
this  paper. 
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Since  a  complete  characterization  of  the  ataoaphere  was  not  available  becauee  of  limitations  in  instru- 
aentation,  it  was  necessary  to  make  soae  improvisations  in  the  coaparison  of  theory  and  experlaent.  For  the 
White  Sands  aeasureaents ,  the  log  amplitude  variance  was  determined  froa  the  chart  recorder  tracings  for  both 
the  94  and  140  GHz  propagation  links.  The  values  of  this  parameter  were  then  substituted  into  the  equation 
for  log  amplitude  variance  derived  by  Hill,  et  al.  to  deduce  the  values  of  the  absolute  hixtldity  structure 
paraaeter  Cq2  which  according  to  theory  would  have  caused  the  millimeter  wave  effects  observed.  This  para- 
aeter  (along  with  C tq!  is  the  only  factor  in  the  equation  not  aeasured  at  White  Sands.  An  assumption  was 
aade  about  the  relationship  between  Cq  ,  CT2,  and  Cfo  that  provides  a  value  for  the  cross-correlatuon  as  will 
be  explained  later.  The  result  of  the  calculation  is  that  plausible  values  for  CQ2  are  obtained.  For  the 
Burlington  aeasureaents,  values  of  the  log  amplitude  variance  were  determined  irom  chart  recorder  traces 
for  four  different  frequencies  .  and  two  different  ranges  L.  For  the  longer  range,  -  2  follows  the  the  theo¬ 
retical  7/6  variation  very  well,  but  the  agreeaent  for  short  path  is  not  as  good.  The  data  follow  the  L1126 
range  dependence  fairly  closely.  Both  sets  of  aeasureaents  will  be  discussed  in  more  detail  in  the  following 
sections. 


Results  of  Theory 

Hill,  Clifford,  and  Lawrence3  have  derived  an  expression  for  the  log  amplitude  variance  of  the  fluctua¬ 
tions  of  electromagnetic  radiation  propagating  through  the  turbulent  atmosphere  which  holds  for  microwave 
through  optical  frequencies  and  maintains  the  dependence  on  Cq2.  CT2,  and  CTq.  This  equation  comprises 
separate  contributions  from  the  real  and  imaginary  parts  and  from  the  cross-correlation  of  these  parts,  as 
follows : 
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1  IK 
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where  K,  1,  and  IK  refer  to  the  real 
part  of  -  2  is  given  by 
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where  k  is  wavenumber,  <T>  is  the  mean  value  of  temperature,  <Q>  is  the  mean  value  of  absolute  humidity,  and 
l*o  1®  tl)e  outer  scale  dimension.  The  dimensionless  parameters  Aj  and  A^  are  given  by 
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where  N  is  the  refractlvlty  of  the  atmosphere.  A  third  equation  involving  A^  where  P  is  atmospheric  pressure 
must  also  be  considered  but  pressure  fluctuations  are  shown  in  Reference  3  to  have  a  negligible  effect  on 
refractive  index  fluctuations  at  all  frequencies  of  interest.  The  refractlvlty  N  is  given  by: 

N  -  (n-11  x  10b  -  (77.6  *K/mb)  |  +  (1.71  x  1U3  *K/g/m3>^  ,  (5) 


where  P  is  in  millibars,  T  is  in  degrees  Kelvin,  and  g  is  in  g/m3.  This  equation  neglects  terms  which  re¬ 
present  dispersive  contributions  to  the  real  part  of  the  refractive  index  in  the  neighborhood  ot  millimeter 
wave  absorption  lines.  References  3  and  6  show  plots  of  Che  refractlvlty  which  support  the  fact  that  dispersive 
contributions  are  at  most  a  few  percent  in  the  frequency  regions  between  lines,  and  therefore  they  would  not 
have  an  important  contribution  to  turbulence  phenomena.  Combining  tquations  (31,  (41,  and  (51  finally  gives 
for  A_  and  A  : 
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The  contribution* 
because  they  contain 


of  the  parameter*  s*  and 
a  parameter  B^  given 1 by: 


0^  to  the  log  amplitude  variance  are  considered  negllgigig 


(5) 


where  B  la  the  ataoapheric  absorption  coefficient.  In  the  low-frequency  window  regions,  tnis  parameter  it 
considered  negligible,  as  is  borne  out  by  calculations  displayed  graphically  in  utterances  3  and  6. 


Finally,  it  is  necessary  to  devise  tone  sort  of  value  for  the  cross-correlation  tent  in  bquatiun  (z;. 
Hill,  et  al.3  suggest  that  an  order-of-magnltude  estimate  for  is 


assuming  a  correlation  coefficient  of  +  1.  For  calculations  in  this  paper,  the  positive  sign  was  cnosen 
because  data  were  taken  on  sunny  days  when  the  surface  was  a  source  ot  both  humidity  and  heat  -  implying 
that  humidity  and  temperature  fluctuations  would  have  a  strong  positive  correlation. 

Using  the  above  results,  it  is  possible  to  substitute  into  kquation  (2)  and  calculate  values  of  Cq^  based 
on  substituting  values  of  .3  2  obtained  from  chart  recorder  tracings  of  field  measurements.  A  later  section 
treats  the  comparison  of  values  of  C  determined  in  this  way  at  VS  and  140  GHz ,  and  also  compares  tne  ire- 
quency  and  range  dependences  of  field  measurements  to  those  given  by  kquation  (21. 

Description  of  kxperlments 

Two  sets  of  measurements  are  used  for  comparison  with  theory  in  this  paper,  namely  the  propagation  tests 
mentioned  earlier  that  were  conducted  at  White  Sands  Missile  Range,  NM,  and  those  conducted  at  Burlington, 

VT.  For  both  sets  of  tests,  meteorological  data  including  temperature,  humidity,  wind  speed  and  direction, 
and  optical  Cn-  were  generally  available;  but  unfortunately  many  of  the  sensors  for  the  Burlington  measure¬ 
ments  were  damaged  in  a  storm  the  night  before  the  measurements  were  made.  Measurements  of  Cq~  and  Gjq  were 
not  available  for  either  case,  which  emphasizes  a  serious  need  for  making  dedicated  measurements  ot  milli¬ 
meter  wave  turbulence  effects.  Both  sets  of  turbulence  data  were  obtained  incidental  to  other  propagation 
tests.  A  system  designed  for  making  such  dedicated  measurements  is  described  in  Reference  7.  j 

Figure  1  is  a  sketch  of  the  experimental  arrangement  used  for  the  white  bands  tests,  which  were  conducted 
at  94  and  ISO  GHz  over  a  2  km  path  in  the  desert.  The  transmitters  used  for  both  frequencies  were  tree- 
running  klystron  oscillators  with  about  30  mw  power  output.  A  precision  rotary  vane  attenuator  was  used  to 
calibrate  several  levels  of  attenuation  for  the  VS  GHz  system,  but  It  was  necessary  to  calibrate  tne  iso  GHz 
attenuation  by  simply  blocking  or  unblocking  the  transmitter.  The  transmitters  were  loo*  modulated  witn  a  1 
kHz  square  wave  to  facilitate  the  use  of  phase  sensitive  detection.  00th  transmitter  and  receiver  antennas 
were  horn-lens  combinations  with  beamwidths  of  2  degrees  for  each  frequency,  and  measurements  were  made 
simultaneously. 

Superheterodyne  receivers  were  used  for  both  frequencies.  The  94  GHz  receiver  used  a  fundamental  mixer 
with  a  local  oscillator  near  94  GHz,  while  the  Ido  GHz  receiver  used  a  harmonic  mixer  pumped  by  a  7u  ghz 
klystron.  Phase  sensitive  detection  with  lock-in  amplifiers  was  used  for  both  cnannels,  with  the  1  knz 
reference  from  the  transmitter  modulator  transmitted  over  twisted-pair  line  to  the  receiver,  oncer  tnese 
conditions,  it  was  necessary  to  reduce  transmitted  power  to  avoid  saturating  the  receiver.  Receiver  outputs 
were  recorded  on  a  dual  channel  chart  recorder. 

TWo  different  systems  were  used  to  make  the  Burlington  measurements.  The  first  system  used  transmitter 
and  receiver  van*  separated  by  about  bUO  a,  and  operated  at  frequencies  of  94,  Ido,  and  22U  ghz.  The  trans¬ 
mitters  were  based  on  the  pulsed  extended  interaction  oscillator  tube  built  by  Varlan  Canada,  and  the  re¬ 
ceivers  were  simple  video  detectors.  Cassegrain  antennas  with  60  cm  diameters  were  used  for  both  transmittml 
and  receivers.  Unfortunately,  the  data  obtained  with  this  system  did  not  have  sufficient  gain  to  permit 
accurate  determinations  of  the  values  of  3^  from  the  recorder  tracings.  Data  from  the  second  system  were 
therefore  used  for  the  comparison  to  be  discussed  in  the  next  section. 

The  second  system  Included  pulsed  radars  st  33,  93,  140,  and  217  Ghz  which  used  IMPATT  transmitters  and 
local  oscillators.  This  system  used  corner  cube  reflectors  at  ranges  of  132  and  bio  m  to  generate  returns 
and  was  originally  designed  to  measure  absolute  atmospheric  attenuation  at  these  frequencies.  Figure  2  is  a 
schematic  diagram  of  the  second  system. 
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The  measurements  aade  at  White  Sands  Missile  Kange  were  divided  into  a  nun  bar  of  aeparate  ‘events'  in 
which  an  explosive  charge  was  detonated  In  the  propagation  path  and  the  subsequent  deterioration  of  milll- 
aeter  wave  signals  was  observed  at  the  receiver.  Measureaents  of  turbulence  effects  were  necessarily 
confined  to  tlaes  slightly  before  and  after  these  events,  because  the  prlaary  purpose  of  the  experlaents  was 
to  observe  propagation  through  dust.  Unfortunataly ,  there  were  only  five  events  In  which  both  the  94  and  140 
GHz  systeaa  operated  aiaultaneouely  well  enough  to  aeanlngfully  coapare  turbulence  effects. 

2 

In  coa paring  theory  and  experiaants,  measured-values  of  ,  P,  T,  and  q  were  substituted  Into  liquations 
(2),  (3),  (4),  and  (3)  to  deteraine  values  of  Cq  for  both  94  and  140  GHz.  Table  1  gives  values  of  para- 
aeters  used  to  sake  these  calculations.  C_2  was  determined  from  aeaaureaenta  of  optical  scintillation  by 
using  the  relation  T 


2 

where  C  la  the  visible  wavelength  index  of  refraction  structure  paraaeter,  and  P  is  aeaaured  in  allllbars. 
Table  IT  gives  values  of  obtained  froa  chart  recorder  tracings  and  values  of  C  2  determined  using 
Equation  (10),  for  each  event.  Simultaneous  data  at  the  two  frequencies  were  not  Obtained  for  Events  A-l  and 
B-l,  as  indicated  in  the  Table.  Figure  3  compares  Cq  determined  in  this  way  from  94  and  14o  GHz  measure¬ 
ments  for  each  of  the  five  events  for  which  simultaneous  data  were  obtained.  The  solid  curve  has  unity 
slope,  and  all  of  the  points  should,  of  course,  lie  on  this  straight  line  for  ideal  agreement  with  theory. 


It  is  also  possible  to  compare  measurements  of  the  frequency  dependence  of  the  log  amplitude  variance  to 
theory.  The  log  amplitude  variances  calculated  at  the  two  frequencies  are  related  by 


°X2(140) 
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according  to  Equation  (2).  Figure  4  shows  the  comparison  between  values  of  this  parameter  measured  at  the 
two  frequencies  for  the  five  events.  Again,  the  points  should  lie  on  the  straight  line  given  by  Equation 
(11). 


The  variance  of  angle-of-arrlval  fluctuations  for  a  spherical  wave  is  given  by 
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where  r  is  antenna  spacing  and  a  is  measured  in  radians.  Although  angle-of-arrlval  fluctuations  were  not 
measured  at  White  Sands,  it  is  of  interest  to  determine  their  magnitude  using  values  of  C  2  calculated  from 
Equatlon2(2).  In  making  this  determination,  note  that  the  factor  in  parentheses  in  this  equation  involving 
C_  ,  Cq  ,  and  C^q  is  just  the  millimeter  wave  C  - .  Figure  5  is  a  plot  of  peak-to-peak  (Or.)  angle  of  arrival 
fluctuations  versus  C  -  using  values  of  this  parameter  derived  at  94  GHz.  Since  ?ft2  in  millimeter  wave 
atmospheric  windows  is  approximately  independent  of  frequency,  by  Equation  (12)  neglecting  dispersion,  the, 
results  in  Figure  5  may  be  taken  as  a  representative  prediction  for  millimeter  wave  windows.  Values  of  C  - 
are  noted  on  the  curve,  although  this  curve  is  not  intended  to  compare  theory  and  experiment.  It  is  simpjy 
Intended  to  show  expected  values  of  angle-of-arrlval  for  the  observed  range  of  Cn*.  An  antenna  spacing  of 
r*J  m  was  used  for  these  calculations. 

2 

That  the  values  obtained  for  C  are  reasonable  can  be  seen  from  the  following  plausibility  arguments.  It 
water  vapor  were  thoroughly  mlxed"so  that  its  density  fluctuations  were  driven  solely  by  temperature  fluctua¬ 
tions,  the  following  relation  would  hold 
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However,  In  the  ceee  of  very  strong  humidity  fluctuations2. 
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Table  II  gives  calculated  values  of  Cn/C  for  each  of  the  events  treated, 
if  this  ratio,  which  Is  0.55  g2/m to  the  average  value  of  <Q>  /<T>Z  * 
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which  Is  a  reasonable  result.  Indicating  that  the  values  of  obtained  in  this  way  are  plausible. 

Finally,  the  frequency  and  range  dependences  of  o^2  measured  at  Burlington  were  compared  to  experiment  at 
35,  95,  140,  and  217  GHz.  Values  of  o^2  obtained  from  the  chart  recorder  t»glngs  are  plotted  as  a  function 
of  frequency  for  the  610  m  path  In  Figure  6.  The  best-fit  curve  having  a  v  depe^egee  was  drawn  through 
these  points  (virtually  a  straight  line).  A  change  In  0X  corresponding  to  (L./L.)1  was  then  scaled 

light  line  was  drawn  parallel  to  the  first!  —  -*  ~  - 


downward  on  the  graph,  and  a  second  straight 


Values 


for  a 


range  of  152  a  were  plotted  on  this  graph.  As  seen  in  the  figure,  the  agreement  with  the  v  dependence  is 
very  good  for  the  longer  path,  but  not  as  good  for  the  shorter.  It  is  possible  that  near-fleld  effects 
contributed  to  errors  in  the  latter  data.  Also,  these  data  were  unaccountably  noisier  than  the  long  path 
data.  However,  Figure  6  shows  that  the  agreement  between  theory  and  experiment  for  the  frequency  and  range 
dependencies  of  is  generally  good. 


Future  Measurements 


A  dedicated  millimeter  wave  turbulence  measuring  facility  with  several  appropriate  meteorological  sensors 
Is  being  assembled  by  Georgia  Tech  and  N0AA  with  support  from  the  I'.S.  Army  Research  Office.  This  system 
will  consist  of  bistatic  transmitter  and  receiver  systems  housed  in  trucks  and  will  probe  all  of  the 
atmospheric  windows  and  absorption  lines  in  the  range  118  CHz  to  340  CHz. 

The  transmitter  will  be  a  phase-locked  klystron  oscillator  with  an  offset  paraboloid  antenna  mounted  in  a 
step-van.  Four  receivers  with  separations  varying  from  1.4  to  10  m  and  pumped  by  a  single  local  oscillator 
will  be  used  to  measure  turbulence  induced  phase  and  amplitude  fluctuations,  from  which  the  atmospheric 
mutual  coherence  function  can  be  derived.  These  results  will  also  give  the  log  amplitude  variance  and  the 
magnitudes  of  phase  fluctuations  from  which  the  angle-of-arrlval  of  the  transmitted  wavefront  can  be 
determined.  The  receivers  will  be  mounted  in  a  semi-trailer  van  and  will  use  horn-lens  antennas.  Power  from 
the  single  phase-locked  klystron  local  oscillator  will  be  distributed  to  the  four  receiver  mixers  by  a  system 
of  optical  beam  splitters  and  beam  waveguides.  The  propagation  path  is  expected  to  be  about  2  km  in  length. 

A  complete  set  of  meteorological  instrumentation  will  provide  atmospheric  data  coincident  with  the 
measured  phase  and  amplitude  fluctuations.  Measurements  of  C_  ,  C  ,  C„,  temperature,  humidity,  wind  speed 
and  direction,  and  aerosol  parameters  will  be  made.  Data  acquisition  and  some  limited  preprocessing  will  be 
accomplished  by  a  microprocessor-based  computer  system,  which  will  simultaneously  store  atmospheric 
parameters  and  turbulence  data  on  tape.  Provision  is  also  being  made  to  record  measured  results  on  chart 
recorders  and  CRT  terminals  to  continuously  monitor  the  data  acquisition  process. 

This  facility  is  expected  to  be  operational  in  the  fall  of  1982  with  frequency  coverage  from  118-174  GHz, 
with  the  other  frequencies  being  added  later.  Reference  7  gives  a  detailed  description  of  this  measurement 
system. 


Cone luslons 


This  paper  compares  theory  to  experiment  bated  on  two  sets  of  millimeter  wave  propagation  measurements 
obtained  with  limited  meteorological  data.  Agreement  in  all  cases  is  generally  good,  but  it  is  again 
emphasized  that  dedicated  millimeter  wave  turbulence  measurements  must  be  made  before  realistic  conclusions 
can  be  drawn.  At  the  present  time,  only  scattered  data  are  available.  The  millimeter  wave  turbulence 
measuring  facility  discussed  briefly  in  the  preceding  section  is  designed  to  meet  the  need  for  more  detailed, 
dedicated  measurements. 
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Due  to  the  nature  of  the  measurements  discussed  In  thl*  paper,  some  care  should  be  taken  In  relating  then 
to  actual  altuations  or  in  ualng  then  to  dealgn  millimeter  wave  ayateaa.  Although  the  agreeaenta  achieved 
have  been  generally  good,  they  Bay  be  fortuitous.  It  does  appear,  however,  that  the  magnitudes  of  the 
effects  are  great  enough  so  that  turbulence  ahould  be  considered  in  the  design  of  ailllaeter  wave  syateas. 
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Figure  1.  Plan  View  of  Experimental 
Arrangement  Used  for  White  Sands  Tests. 
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Figure  2.  Schematic  of  Burlington  Tests. 


Figure  3. 


Comparison  of  C  ^  Determined  from  94  and 
140  GHz  Measurements. 


Figure  4.  Comparison  of  Determ 
from  94  and  140  GHz  Meas 
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Figure  5 


Calculated  Peak-to-Peak  Angle-of-Arri val 
Fluctuat ions. 


Figure  6.  Comparison  of  the  Frequency 
and  Range  Dependence  o! 
Turbulence  Measurements  to 
Theory . 
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